Periodic Density of States Modulations in Superconducting Bi 2 Sr 2 CaCu208+5 
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In this paper we show, using scanning tunneling spectroscopy (STS), the existence of static striped 
density of electronic states in nearly optimally doped Bi2Sr2CaCu20s+i in zero field. The observed 
modulation is strongest at roughly half the superconducting gap energy and is aligned with the 
Cu-0 bonds, with a periodicity of four lattice constants, exhibiting features characteristic of a 
two-dimensional system of line objects. These features also exhibit asymmetries between the two 
orthogonal directions locally, pointing to a possible broken symmetry state (i. e., stripe phase). We 
further show that the density of states modulation manifests itself as a shift of states from above to 
below the superconducting gap. The fact that a single energy scale (i. e., the gap) appears for both 
superconductivity and stripes suggests that these two effects have the same microscopic origin. 

PACS numbers: PACS numbers: 74.72. Hs, 74.50.+r, 74.25.-q 



I. INTRODUCTION 

One of the most important problems for condensed 
matter physics is the determination of the ground 
states of strongly correlated electron systems. In 
particular, high-temperature superconductors (HTSC) 
appear to have unusual ground states nr^flpfrt^e of 
strong electrgn^jX^jdjaidaas. TheoreticaltrnuaBB and 
experimentafll3'BE3liiJlljllj'0 evidence has been mount- 
ing in support of the possibility that their ground state 
exhibits spin and charge density waves (SDW and CDW), 
which are primarily one-dimensional (i. e., stripes) 
and which either compete with or promote supercon- 
ductivity. Coexistence of CDW or SDW and supercon- 
ductivity Jias^pAdously been reported in the lower T c 
matcrialsQ'PtLOESEJ and in the presence of large mag- 
netic fieldstSO. 

In inelastic neutron scattering experiments on 
La2- x Sr x Cu04 near optimum doping (x=0.163), aiid 
in the presence of a magnetic field, Lake et aliJH 
found strong scattering peaks at the four fc-space 
points: (tt/oo) [(1/2,1/2) ±5(0,1/2)] and (tt/oo) [(1/2,1/2) 
±5(1/2,0)], where S ~ 0.25 and ao is the lattice constant. 
This implies a spin-density of periodicity ~ 8ao, which 
the authors find extends to more than 50aft r -.imich be- 
yond the vortex core. Also, Mitrovic et alS3 reported 
similar periodicity in a high-field NMR imaging exper- 
iment. Using STM, Hoffman et alS3 reported scanning 
tunneling spectroscopy on L^S^CaC^Os+a single crys- 
tals revealing a checkerboard of quasiparticle states with 
four unit-cell periodicity surrounding vortex cores. In 
/c-space this corresponds to Bragg peaks in the local den- 
sity of states at (27r/a )(±l/4,0) and (27r/a )(0,±l/4). 
This structure was found around the energy (~7 meV) 
of a feature previously observed in the spectra of vortices 



in the same field rang The checkerboard structure 

extends over ~r2f)ao, again indicating ordering outside 
the vortex corestll. 



The absence of data showing stripes in the higher T c 
materials in zero field has lent support to the idea that 
stripes are competing with high-temperature supercon- 
ductivity. However, we showedllS, using scanning tun- 
neling spectroscopy (STS), a first indication of the ex- 
istence of static striped density of electronic states in 
nearly optimally doped Bi2Sr2CaCri208+5 in zero field. 
This modulation is aligned with the Cu-0 bonds, with a 
periodicity of four lattice constants, and exhibits features 
characteristic of a two-dimensional system of line objects. 
We further showed that the density of states modulation 
manifests itself as a shift of states from above to below the 
superconducting gap. The fact that a single energy scale 
(the gap) appears for both superconductivity and stripes 
suggests that these two effects have the same origin. In 
the present paper we expand on our initial results, study- 
ing further the energy dependence of the modulation and 
its interrelation with quasiparticle scattering. Our main 
result in the present paper is that the four-period peaks 
in the Fourier transform of the data are present at all 
energies, including very low energies. At the same time, 
peaks at other fc-vectors appear at higher energies, sug- 
gesting the contribution of quasiparticle-like scattering 
from impurities, as suggested by Wang and Leeta. The 
two effects are separable, suggesting that a simple model 
of two-dimensional band electrons cannot explain all the 
observed physics, corrobajjatHig the photoemission find- 
ings on the same systemcSIO'ESl. 
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II. EXPERIMENTAL 

STS allows one to measure, on an atomic scale, the 
electronic density of states. This measurement of the 
local density of states (LDOS) makes STS a powerful 
tool for investigating correlated electron systems. In 
this study we used r^S^CaC^Os+a because it has 
a high T c (^90 K), and since it cleaves easily, yielding 
large, stable, atomically flat surfaces. The single crys- 
tals of Bi2Sr2CaCu208+5 used in this study, grown by a 
floating-zone method, were annealed to be slightly over- 
doped, yielding a superconducting transition tempera- 
ture of 86 K. The samples were cleaved (between the BiO 
planes) at room temperature in a vacuum of better than 
1 x 10~ 9 torr, then transferred in less than one minute 
to the low-temperature STM, where cryopumping yields 
orders of magnitude better vacuum. Data were taken at 
8 K with an iridium tip at a sample bias of -200 mV 
and a setpoint current of -100 pA. The choice of setpoint 
establishes the relatively arbitrary normalization condi- 
tion for the raw differential conductance (dl/dV), which 
is proportional to the local density of states. The data 
sets used in this paper consist of 64 x 64 pixel images 
with dl/dV spectra taken at each point. In addition, we 
took 256 x 256 pixel topographic scans over the same ar- 
eas. Fig. la shows the topography of a typical area: the 
surface exhibits clear atomic resolution. In particular, 
the superstructure in the BiO planco with average pe- 
riodicity ^27 A(as well as the location of the individual 
Bi atoms) is evident in all scans and provides a-xefer.- 
ence direction for our study. As reported earlicrcfloEj, 
all samples also exhibit strong inhomogeneity in the local 
spectroscopy. This manifests itself here as large variation 
in the superconducting gap size (A), as measured by the 
location of the peak in the local density of states (Fig. 
lb). Overall the spectra evolve from small superconduct- 
ing gaps (A ~ 30 mV) to very large gaps (A > 60 mV) 
with mucli—ijpdptjcd coherence peaks reminiscent of the 
pseudogapEftEZlc3. 



III. RESULTS AND ANALYSIS 

A. Fourier Analysis 

Fig. 2a shows a map of the LDOS for 15 mV sam- 
ple bias. While there are some features in the real- 
spacej-image, they are more readily observed in Fourier 
spaceE3. Therefore, we will focus for the moment on Fig. 
2b, which shows the power spectrum, namely the am- 
plitude squared of the Fourier transform, of the data 
in Fig. 2a. Fig. 2b exhibits a large amount of weight 
at small wavevectors, coming mostly from the random 
contributions of the inhomogeneous gap. In addition, 
there are two bright peaks appearing along the k x di- 
rection which originate from the superstructure. Finally, 
Fig. 2b clearly shows four peaks, oriented 45 degrees to 
the superstructure and at points corresponding to, within 



experimental uncertainty, periodicities of (±4ao, 0) and 
(0, ±4a ), oriented along the copper-oxygen bonds. The 
real-space data (Fig. 2a) show a two-dimensional pattern 
in the LDOS that is oriented along the diagonals; how- 
ever, as is evident from Fig. 2b, contributions from fea- 
tures at other wavelengths obscure this modulation. The 
period of the modulation is shown clearly in Fig. 3, the 
autocorrelation function of Fig. 2a. This figure has been 
rotated approximately 45 degrees clockwise from Fig. 2a 
to place the copper-oxygen-copper directions horizontally 
and vertically. The maxima are separated by ~ 4ao . The 
apparent rhombic distortion of the central nine maxima 
comes from the near coincidence of the upper left and 
low right corners with the superstructure period, which 
moves these maxima to slightly larger distances. 

In order to further elucidate the spatial structure of 
the observed modulation, we can filter from Fig. 2b the 
contributions to the LDOS modulation that are far away 
from the (27r/4o )(±l,0) and (27r/4a )(0,±l) points. In 
general, filtering in reciprocal space is performed by mul- 
tiplying the complex two-dimensional Fourier transform 
F(k x , ky) of the real space image f(x, y) with a real filter 
mask m(k x ,k y ) and subsequently calculating the inverse 
transform. The filter mask must be a symmetric func- 
tion with respect to k=0. In order to avoid artifacts, one 
should take care that there is little weight around the 
edges of the filter, aiid, that the edges of the filter are 
not particularly sharpEJ. In addition, examination of the 
resulting image for a variety of widths can show whether 
artifacts are important since the wavevector of such arti- 
facts will depend on the sizes of the filter. For our data, 
we have chosen a radial Gaussian mask function with ra- 
dius a centered around the four-period points. Obviously 
such a procedure will pick out a preferred wavelength; its 
validity here rests on the observation of relatively large 
weight in these regions in the raw power spectrum, and 
the relatively small weight at the filter edges. This pro- 
cedure is hardly novel: similar Fourier filtering was suc- 
cessfully used to remove the topographic ^background in 
order to enhance the Friedel oscillationsE3 of Be(1010) 
STM imagesEJ, for example. 

Fig. 4 shows the real-space data after such Fourier fil- 
tering with a standard Gaussian filter function using a— 
(27r/15ao). Varying the width of this filter has no qual- 
itative effect on the resulting image, and the size of the 
above filter implies about 30 pixels within the FWHM 
around the center of each Bragg-spot. The modulation 
partially visible in Fig. 2b is now clear: Fig. 4 shows a 
dominant four-period modulation of the density of states 
that is almost checkerboard like, but with defects char- 
acteristic of a two-dimensional four-fold structure. In 
particular we observe dislocations in the form of extra 
half-rows. These defects locally increase or decrease the 
periodicity. From the separation of the defects in this 
map we estimate the correlation length to be ~ 20ao. 

The defects noted above in Fig. 4 can also explain the 
weight distribution of the four-period peak in the Fourier- 
transform. Given the approximately 80 Acorrelation 
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length of the modulation, we would expect a peak width 
of about 27r/20ao. However, because the transform is 
taken over an area that is only about two correlation 
lengths wide, we expect to see contributions at only some 
of the wavevectors in this region. Since any dislocation 
alters the local periodicity of the checkerboard, we see 
a number of peaks at positions near 2-7r/4ao. Assum- 
ing that these reflect a single peak broadened by the fi- 
nite correlation length, we therefore place the peak at 
27r/4a ±27r/40a . 



B. Energy Dependence 

To further understand the nature of the modulation 
near k x . y — 27r/4ao, we studied the energy dependence of 
the effect. Fig. 5 shows the power spectra of the conduc- 
tance over this area for a number of voltages. As above, 
the axes have been rotated to place the Cu-O-Cu direc- 
tions parallel to the axes. For all sub-gap energies, there 
is significant weight in the vicinity of the 27r/4ao points, 
as shown by the circles. However, the small wavevector 
(long wavelength) contributions increase with increasing 
energy, up to the nominal gap size, as we would expect 
from the spatial variation in gap size. This means that 
the 27r/4ao points are most relatively prominent for lower 
energies. 

The four-period peaks are not the only features ob- 
served in the power spectra. Also evident are the su- 
perstructure points at ~ 27r/ao(±0.1, =f0.1). In addi- 
tion, for most energies one can see faint second-order 
superstructure spots. Curiously, there also appear spots 
perpendicular to the superstructure at the points ~ 
27r/ao(±0.1, ±0.1), which do not appear to disperse with 
energy. Finally, while there is a significant weight at 
many wave vectors within the box bounded by the k x , y — 
2-7r/4ao points, there is, with the exception of the second 
order superstructure spots, essentially no weight outside 
of it. 

Fig. 6 shows a number of line cuts along the (tc, 0) and 
(0, 7t) directions. Each curve shows the magnitude of the 
Fourier transform along the Cu-O-Cu direction for a dif- 
ferent energy. The curves have been computed by Gaus- 
sian averaging in the vicinity of the desired point with a 
FWHM of one pixel width to remove pixel size artifacts 
and account for the (small) uncertainty in knowing the 
Cu-O-Cu direction. The peaks around k x ^ y = 2ir/Aao 
show up at all these sub-gap energies in both biases 
(although only the positive bias is shown here for clar- 
ity). While there appears to be a small change, from 
0.25 x (27r/ao) at low energy to about 0.23 x (27r/ao) 
around 30 mV, no significant dispersion is found for these 
peaks. Whether this small shift comes from the appear- 
ance of an additional dispersive or non-dispersive feature, 
an increase in the low- wavevector background, or a small 
dispersion in these k x , y ~ 27r/4ao points should be re- 
solved by additional data. 

In a recent theoretical paper, Wang and LeeEj ex- 



plored the possibility that band-structure quasiparticles 
will exhibit quantum interference due to the presence of 
point defects in the system. This would manifest itself 
as strongly dispersing features in the Fourier transform 
of the LDOS as a function of energy. Clearly, such an 
effect will be strong at high energies (bias voltages) and 
will disappear close to zero bias due to the lack of phase 
space for the possible k's. Indeed, Wang and Lee find 
essentially no peaks at any k for energies close to zero 
bias in their theoretical study. Since we see a feature at 
0.25 x (27r/ao) for V = mV and it shows no dispersion 
up to 20 mV, it is clear that we observe something other 
than the dispersion predicted by Wang and Lee. We will 
come back to this point in the next section. 



C. Phase Information 

In addition to the information extracted above from 
the magnitude of the Fourier transform, there is also in- 
formation contained in the phase for the variety of en- 
ergies sampled. We now use this to further understand 
the four-period modulation observed and contrast it with 
other features discussed above (Figs. 5,6). Fig. 7 shows 
that the modulation in the LDOS exhibits a strong en- 
ergy dependence. Here we show the value of the Fourier 
transform at the two independent peaks from Fig. 2b as 
a function of the bias voltage. The error bars are deter- 
mined by modeling the variations in dI/dV(x,y) at each 
energy by uncorrelated noise of the same amplitude. The 
arbitrary phase (dependent on the point about which the 
Fourier transform is computed) has been set to maxi- 
mize the real part, plotted in red, which changes sign at 
about ±40 mV. The variation in the imaginary part of 
the signal is considerably smaller, consistent with zero 
within the uncertainty. This is important, since it shows 
that the location in real-space of the density of states 
(DOS) modulation is the same at all measured energies, 
although at about ±40 mV the positions of the maxima 
and the minima switch. 

The constancy of the phase can be illustrated in an- 
other manner, reminiscent of the analyses of neutron 
scattering results: Fig. 8 is similar to Fig. 5, but in- 
tegrated in energy. Each panel shows the power spec- 
trum of the current at that energy, which is equal to the 
integral from zero to V of dl/dV. It is clear that with 
increasing energy, the four-period modulation remains 
relatively strong, while many other features fade away. 
The features that remain at many energies are appar- 
ently summed with the same wavevector and phase, i. e., 
they are a non-dispersive and pinned. On the other hand, 
many of the other features (internal to the four spots in 
fc-space) are reduced in intensity, indicating more varia- 
tion in either phase or fc-space location with energy. To 
focus on the peaks at k x>y = 27r/4ao, we show in Fig. 9 
the contrast between the raw dl/dV and the integrated 
dl/dV for line scans along the (0,0) to [tt, 0) direction 
as a function of wave-vector and energy. The fact that 
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the peak at k x = 27r/4ao is enhanced as a function of en- 
ergy for the integrated cupwes reflects the non-dispersing 
nature of this modulationElJ. 

Furthermore, the energy dependence of the Fourier 
transform (Fig. 7) shows that relative spectral weight 
is shifted to sub-gap energies (with a peak at ±25 mV) 
from intermediate energies (with weight between 50 and 
150 mV). The comparison of the above result with the 
energy of the superconducting gap (A ~ 40 mV) provides 
strong circumstantial evidence that the striped quasi- 
particle density and superconductivity are intimately 
connected^. Since the integrated density of states from 
zero to infinity is the total charge, the shift of weight from 
intermediate to low energies coupled with the small mag- 
nitude of the Fourier components (< 5% of the LDOS) in- 
dicates that the total charge participating in this density 
of states modulation is very small, as is expected : 
erant systems with strong Coulomb interactionsESI 
On the other hand, it is important to note that energy 
independent spatial variations in the DOS would not ap- 
pear in this analysis because the normalization is set by 
maintaining constant current (which is proportional to 
the integrated DOS from zero to the setpoint voltage). 
We are sensitive to features in the shape of the LDOS, 
not the overall magnitude. Energy independent changes 
would however appear in the so-called topographic sig- 
nal (Fig. la), which actually maps a contour of constant 
integrated density of states. The power spectrum of the 
topograph reveals no such peaks above the noise; how- 
ever, this is inconclusive because of the small amount of 
total charge exe<pd©d to participate in the modulation 
for this systemKJt 



D. Anisotropy 

To examine whether the LDOS modulation in our ex- 
periments is anisotropic, we performed Fourier analysis 
on many smaller sub-regions of a slightly larger area (Fig. 
10). In Fig. 10c and lOd we show maps of the spatial 
variation of the amplitude of the two perpendicular pairs 
of Bragg peaks. The color at each point is determined 
by taking the magnitude of the Fourier transform at the 
four period point for the 32 x 32 pixel dI/dV(15mV) map 
centered at this point. In addition to the small scale 
variations which come from noise and finite size effects, 
it is clear that each varies on the scale of roughly 100 
A. In addition, these large scale features exhibit signif- 
icant anticorrelations between the two maps: there are 
several regions in space where one map shows a local 
maximum while the other shows a local minimum. This 
indicates that locally there is a broken symmetry, and 
that the underlying order is not two-dimensional, but 
one-dimensional. The two-dimensional nature of the pat- 
terns observed appears to come from the interpenetration 
of these perpendicularly oriented stripes. 



IV. DISCUSSION 

"Stripe order" usually describes unidirectional density 
wave states, which can involve unidirectional charge mod- 
ulations ( "charge stripes" ) or coexisting charge and spin 
density order ("spin stripes"). Such a state in HTSC, 
with spatial segregationjpf charge, was ficst suggested by 
Zaanen and GunnarsonEI and by Schulzo. Stripe order 
is expected to occur via a quantum phase transition (i. 
e. at T — 0), where the broken symmetry state is one 
that breaks rotational symmetry and translation symme- 
try perpendicular to the stripe direction (in a crystal we 
mean breaking the crystal symmetry group.) For HTSC, 
in the regime where superconductivity is present, Kivel- 
son and Emeryo argued that while no long range stripe 
order will be present for T > 0, the proximity of the sys- 
tem to the stripe quantum critical point implies strong 
stripe fluctuations. In this stripe- liquid phase stripes 
are dynamical objects and thus cannot be detected with 
static probes such as STM. However, even a small amount 
of quenched disorder will pin the stripes, hence re 
ing the nature of the nearby ordered stripe phaseE 
VortejCj-induced pinning of stripes was used by Hoffman 
et oZ£a to explain the observed density modulation near 
vortex cores. Random point defects, present in all HTSC 
samples, could play a similar role, thus explaining the 
origin of the stripe patterns we observe. 

While the analysis presented above concentrated on 
the robust feature at k- vectors 27r/4ao along the Cu-0 
bonds, it is evident that our data show more than just 
that modulation. Features, which may bedispersing sim- 
ilar to the suggestion of Wang and LeeEj, appear in the 
data. Thus, while it is clear that the simple picture of 
Wang and Lee cannot explain our data, the extra features 
suggest that the stripe picture alone is not comprehen- 
sive. Explaining the electronic structure data using a 
simple band structure approach was aheadyjhnw.p to be 
inadequate by recent ARPES studiesBBHOHEl. In 
the normal state (i. e. at temperatures above T c ), the 
main problem arises in analyzing the single hole spec- 
tral function as a function of frequency and wavevector. 
While for a fixed k the spectral function showed only 
broad features as a function of lo, the structure of the 
spectral function became sharp when measured at a fixed 
frequency u, and as a function of k. This effect points 
to the absence of stable excitations with the quantum 
numbers of the electron, hence, a non-Fermi liquid situ- 
ation exists in the normal state of the cuprates. In the 
ordered state (the superconducting state) quasiparticles 
seem to become well definedc3. This fact was taken at 
face value by Wang and Lee to suggest that in the super- 
conducting state they should be able to use the simple 
band structure approach. 

However, the same ARPES data clearly indicate that 
the superconducting state is anomalous as well, with a 
coherence-like peak that tracks the superfluid density and 
the hole doping instead of -the gapc3, and an anomalous 
nodal quasiparticle lifetimes. Here it was found that the 
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inverse lifetime of the excitations is proportional to tem- 
perature, independent of binding energy, for low energies, 
and is proportional to energy, independent of tempera- 
ture for large binding energies. This "marginal Fermi 
liquid" behaviorE-J is unaffected by the occurrence of su- 
perconductivity as the temperature is lowered. Never- 
theless oscillations that are of a single particle nature 
will occur even when there are no well defined quasi- 
particles, so long as there are some elementary excita- 
tion in the system with a well defined dispersion rela- 
tion. If we take the underlying band structure measured 
by ARPESpfjau. Bi 2 Sr 2 CaCu208+5 at similar doping as 
a guidelinetiiE3, we find discrepancies with the quasi- 
particie scattering interference picture in two regimes. 
First, as we noted above, we observe Fourier peaks at 
0.25 x (27r/ao) at and near zero bias. It is easy to see 
that the quasiparticie interference picture cannot pro- 
duce such peaks at low energieaLjCJ. Assuming quasipar- 
ticles with fc-dependent energy: E^ = •v/A? + et, where 
Afc = ^[cos(k x CLo) ~ cos(kyao)] is the d-wave supercon- 
ducting gap and is the normal-state quasiparticie dis- 
persion (assuming well-defined quasiparticles exist) mea- 
sured from the Fermi surface. We assume scattering on 
the Fermi surface from state k to k' with q = k' — k, such 
that q = (27r/4ao,0). In that case k = (— ir/Aao, k y ), 
and k' = (tt /Aa a ,.k y ). If we extract k y directly from the 
ARP^p^ataBB, we find k y 0.6(7r/a ). For A - 30 
meVE3c3'cZI, which is also the minimum gap found on 
this type of samples using STMuH, this gives an esti- 
mate of the lowest energy for quasiparticie scattering at 
this wavevector of approximately Ao/2 ~ 15 meV. The 
correct energy cutoff may be somewhat different depend- 
ing on the exact details of the band-structure, but this 
analysis certainly excludes the energies around zero bias. 
Second, as noted above, the assumption of quasiparticles 
with well defined dispersion should cease to be justified 
when the energy is ~ fcgT c ~ 15 meV since at these en- 
ergies we expect the behavior to be more normal-state 
likec3. This is in particular true if the system is closa-to 
a quantum critical point as was shown by Valla et a/.a. 

Our results however are consistent with an underlying 
stripe- structure, as are, in a similar way, the ARPES 
datallj. If we adopt a quasi one-dimensional approach 
(i. e., stripe structure) to describe the normal state, it 
is clear that since continuous global symmetries cannot 
be spontaneously broken in one dimension, even at T = 
0, the systemJias to exhibit interchain coupling at the 
onset of ordeiCJ. In that context the system will retain 
some of its one dimensional nature, and thus one should 
observe in STM studies of point impurities in high-T c 
the signature of both Friedel-like oscillations and a fixed 
CDW at the stripe wavevector. 

Modulations which reflect the spectrum of elementary 
excitations can however be distinguished from those re- 
lated to incipient order. In particular, the peak in the 
Fourier transform of the LDOS produced by incipient or- 
der tends to be phase coherent, in agreement with Fig. 7. 
At the same time, features that are not due to incipient 



order will either have a random phase (if they originate 
from noise) , or a phase that is strongly energy dependent 
(if they originate from dispersing features). Thus, as sug- 
gested in the previous section, in order to focus on those 
aspects of the STM spectrum that reflect incipient order, 
one should integrate over a small but finite range of ener- 
gies about E = 0; non-dispersive, phase coherent features 
of the spectrum will thus be accentuated, while disper- 
sive, phase incoherent features will tend to canceEJ as 
was demonstrated in Figs. 8 and 9. Moreover, the energy 
dependence presented in Fig. 7 was obtained ttepsetically 
for models that include striped-charge orderESBa. All 
these models rely on pinning of dynamic spin and charge 
density waves fluctuations in a d-wave superconductor 
by local imperfections which preserve spin-rotation in- 
variance, such as impurities or vortex cores. In partic- 
ular, Podolsky et al. derived an explicit expression for 
the Fourier component of the local density of states at 
the ordering wavevector, and demonstrated that various 
patterns of translational symmetry breaking exhibit dif- 
ferent energyj-dependence. Comparing to our data pre- 
sented earlieita, which is reproduced here in Fig. 7, they 
argued that the data are consistent with weak transla- 
tional symmetry breaking that includes charge density 
waves and direct dimerization. Indeed Fig. 4 in their 
paper reproduces remarkably well our results. 

Finally we note that in addition to this static striped 
density of states structure that spans the entire sur- 
face, gap structure inhomogeneities are observed in all 
our samples. The movement of LDOS to the middle of 
the gap, around a particular energy, and its origin from 
states just above the gap suggests that the density mod- 
ulation and superconductivity are inseparable. The inho- 
mogeneities in the superconducting gap structure should 
therefore also be reflected in the striped density of states. 
Indeed, the striping does show local variation as is shown 
in Fig. 9, together with topological defects in the striped 
structure visible in Fig. 4EJ. 



V. CONCLUSION 

In this paper we studied the LDOS spectra of slightly 
overdoped BSCCO single crystals. Based on our data 
we conclude that the studied samples exhibit two types 
of modulated features in the energy dependent LDOS. 
First, there are peaks at wavevectors corresponding to 
four-period modulation of the LDOS along the Cu-0 
bond direction, which we take as evidence for stripe 
order. Second, there are features, primarily at shorter 
wavevectors, which may to be energy dependent. These 
may be attributed to a combination pf elementary ex- 
citation (quasiparticle-like) dispersiontj and the strong 
inhomogeneities i^thft gap structure previously observed 
by several groupsLaOES. Finally, it is remarkable that in 
this single set of experiments we see locally the d-wave 
superconducting gap, the pseudogap, and metallic 
stripes, suggesting that all reflect aspects of the same 
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physics. 

While this paper was being completed, Hoffman 
et alE3 published higher resolution but qualitative^ 
similar data to that first published by Howald et al$3 
and to the data presented in this paper. However, in 
their paper they focus on the dispersive features of the 
data, overlooking the existence of peaks at the stripe 
wavevector which appear for all their reported energies. 
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FIG. 1: Inhomogeneity in the electronic structure of a slightly 
overdoped Bi2Sr2CaCu208+* single crystal, a) Topography 
(160 Ax 160 A) of the cleaved BiO surface, b) Gap size over 
the same area. 
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FIG. 2: Periodic spatial variation in the quasiparticle density 
of states, a) Map of dI/dV(15 mV). b) Power spectrum of 
the Fourier transform of dI/dV(15 mV). 
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FIG. 3: Autocorrelation function of Fig. 2a showing the pe- 
riod of the primary modulation (units are ao). Axes have 
been rotated from Fig. 2a to make the Cu-O-Cu directions 
horizontal and vertical. 
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FIG. 4: Fourier filtered real-space map of dI/dV(15 mV) (Fig. 
2a), using a sum of Gaussian fc-space filters centered around 
k x>y — ±27r/4ao. 
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FIG. 5: Power spectra of dl/dV as a function of energy. Energies are 0, 3, 7, and 10 mV (top row); 13, 16, 19, and 22 mV 
(middle row); 25, 28, 32, and 35 mV (bottom row). White circles show the width (a) and position of the Gaussian filter used 
in Fig. 4. Axes have been rotated from Fig. 2. 
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FIG. 6: Line scans along the (0,7r) and (tt,0) directions, re- 
spectively, through each of the two independent Fourier trans- 
form peaks as a function of energy (\dl /dV(0, k y , E)\ and 
\dl /dV(k x , 0, E)\). Scans are shifted vertically for clarity. 




FIG. 7: Fourier transform at k=(27r/4a )(±l,0) (top) and 
k=(2-7r/4ao)(0,±l) (bottom), the location of the peaks in Fig. 
2b, as a function of sample bias. The red and blue traces 
correspond to the real and imaginary parts, respectively. 
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FIG. 8: Power spectra of I (integral of dl/dV) as a function of energy. Energies used are the same as Fig. 
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FIG. 9: Line scans as a function of k x along the (0, 0) to (71", 0) 
direction, and as a function of energy (color scale). Top panel 
shows the LDOS (dl/dV), and bottom panel the integrated 
LDOS up to the given energy. 
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FIG. 10: Spatial variation in the amplitude of the periodic 
modulation over a larger area, a) Constant current topograph 
(130 Axl30 A), b) Power spectrum of dI/dV(15 mV) over the 
entire 260 Ax 260 Aarea showing the peaks used in c) and d). 
c) Magnitude of the local Fourier transform of dI/dV(15 mV) 
at one of the 2n/4ao points for a number of 32 x 32 pixel 
regions, d) The same map as c), for the other peak. 



